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1.  BACKGROUND  AND  PROJECT  OBJECTIVES 


During  a  two-week  period  in  July  1990,  the  United  States  Air  Force  Phillips  Laboratory 
(GP)  conducted  an  intensive  field  experiment  using  a  variety  of  optical  sensors  at  the  Brunswick, 
Maine,  Naval  Air  Station.  This  trial  was  named  FLAPIR  (an  acronym  for  R  TR*  and  lidar  atmo¬ 
spheric  propagation  in  the  infrared).  The  objective  of  the  FLAPIR  program  is  to  determine 
whether  uiformation  obtained  from  lidar  measurements  can  be  used  to  predict  the  performance  of 
a  forward-looking  infrared  (FLIR)  electrooptical  system  during  adverse  weather  conditions.  In 
support  of  Phillips  Laboratory,  SRI  International  operated  its  four-wavelength  lidar  system  dur¬ 
ing  the  FLAPIR  field  tests.  TTie  primary  role  of  SW’s  participation  was  to  operate  the  lidar  sys¬ 
tem  and  to  provide  Phillips  Laboratory  with  the  lidar  backscatter  data  obtained  during  the  trial. 
SRI  was  also  to  assist  Phillips  Laboratory  in  the  analysis  of  those  data  by  calculating  path 
transmission  values  for  comparison  with  the  dual-wavelength  infrared  (DUWIR)  camera 
measurements. 


•An  acronym  for  forward-looking  infrared. 
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2.  SRI  FOUR-WAVELENGTH  LIDAR  SYSTEM 


Under  funding  from  the  Army  Research  Office,  SRI  developed  a  four-wavelength  lidar  sys¬ 
tem  to  investigate  multiple-wavelength  lidar  techniques  for  quantitative  aerosol  property  mea¬ 
surements  from  remote  distances.  The  system  is  also  effective  for  evaluating  propagation  effects 
at  wavelengths  commonly  used  by  DoD  laser  systems.  The  SRI  system,  which  is  installed  in  a 
20- ft  van,  transmits  laser  pulses  at  four  wavelengths:  0.53  and  1.06  pm  using  an  Nd  laser, 

3.8  pm  using  a  DF  laser,  and  10.6-p.m  using  a  CO2  laser.  An  optical  arrangement  provides  for 
transmitting  the  pulses  along  the  same  horizontal  path,  coaxial  with  a  common  12-inch-diameter 
Newtonian  receiver.  Wavelength-separation  optics  direct  received  energy  to  separate  detectors 
for  0.53-pm,  1.06-pm,  and  3.8/1 0.6-pm  wavelength  energy.  A  receiver  multiplex  circuit  routes 
the  detected  0.53-pm  backscatter  signature  to  a  transient  recorder  on  the  first  Nd  laser  firing,  and 
the  1.06-pm  backscatter  si^ature  on  the  second  Nd  laser  firing.  A  double-pulse  mode  provides 
both  firings  within  a  50-ps  period.  The  multiplex  unit  then  switches  to  the  3.8/10.6-pm  receiver 
and  the  DF  laser  is  pulsed  about  15  p.s  later  followed  by  the  CO2  laser  firing  at  another  15-ps 
interval. 

The  transient  recorder  is  controlled  by  an  interrupt  circuit  that  terminates  the  digitization  of 
each  firing.  This  arrangement  provides  for  storage  of  a  four- wavelength  backscatter  signature 
within  the  transient  recorder  within  a  80-p.s  interval.  The  logarithmic  backscatter  values  are  then 
transferred  to  an  LSI-1 1/24  microcomputer  and  written  on  nine-track  magnetic  tape,  together 
with  date  and  time  information. 
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3.  FIELD  TEST  AND  DATA  COLLECTION 


During  the  FLAPIR  field  trial,  the  SRI  four- wavelength  lidar  van  was  posidoned  alongside 
other  instrumentation  vans  at  the  Brunswick  Naval  Air  Station  site.  A  passive  target  was  placed 
at  a  distance  of  928  m  from  the  lidar  van  to  provide  for  lidar  sighting,  calibration,  eye-safety,  and 
target  returns  useful  for  data  analysis.  The  Udar  was  operated  during  assigned  time  periods 
between  13  July  and  27  July  1990.  The  data-acquisition  periods  on  13  July  served  as  a  system 
shakedown.  Data  at  each  wavelength  were  digitized  at  a  sampling  interval  of  20  ns  with  a  preci¬ 
sion  of  eight  bits,  and  all  data  were  recorded  on  nine-track  magnetic  tape.  A  real-time  color 
graphics  display  was  used  to  evaluate  data  quality  and  prevailing  atmospheric  scattering  proper¬ 
ties  at  each  wavelength.  No  major  equipment  problems  occurred  during  the  period  16  July 
through  27  July,  although  several  computer  failures  resulted  in  short  data  gaps. 
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4.  DATA  REDUCTION  AND  ANALYSIS 


A.  DATA  STORAGE  AND  TRANSMITTAL 

At  the  time  of  data  acquisition,  the  lidar  backscatter  signatures  at  the  four  wavelengths  were 
recorded  together  with  time  and  date  information  on  nine-track  magnetic  tape  as  a  single  fixed- 
length  record  of  8256  eight-bit  bytes  for  each  record.  The  actual  backscatter  data  were  stored  in 
the  last  8192  bytes  (4096  16-bit  words).  These  data  were  subsequently  copied  to  a  MicroVAX 
computer,  where  they  are  stored  on  optical  disk  as  unformatted  direct  access  files  (in  VaxA^MS 
Files- 1 1  format)  containing  fixed-length  (8256-byte)  records,  with  a  separate  file  for  each  physi¬ 
cal  raw  data  tape.  Each  data  tape  could  hold  up  to  1200  lidar  records;  thus,  the  maximum  num¬ 
ber  of  records  per  file  is  1200,  although  the  actual  number  varies. 

An  optical  disk  containing  the  raw  lidar  data  files  in  VAX/VMS  Files- 1 1  format  has  been 
sent  to  Lt.  Col.  George  G.  Koenig  at  Phillips  Laboratory.  This  disk  contains  all  data  from 
16  July  through  27  July.  Data  collected  on  13  July  were  excluded  because  no  DUWIR  data  were 
taken  on  this  day.  A  simple  FORTRAN  computer  program  for  reading  the  data  files  was 
included  on  the  disk,  and  written  instructions  were  provided.  At  Col.  Koenig’s  request,  a  dupli¬ 
cate  disk  was  sent  to  Dr.  Luc  Bissonnette  of  Defense  Research  Establishment  Valcartier. 

B.  TARGET  IDENTIFICATION  AND  DETECTABILITY 

The  first  data-reduction  task  was  to  examine  each  file  and  to'identify  when  the  passive 
reflective  target  was/was  not  detectable  at  each  wavelength.  This  multiple-step  procedure  was 
described  in  detail  in  the  first  SRI  Quarterly  R&D  Status  Report  (15  October  1990).  Basically, 
the  procedure  required  graphic  display  of  numerous  relative  backscatter  signatures,  and  it  com¬ 
bing  subjective  analysis  with  an  objective  search  routine  to  identify  the  presence  or  absence  of  a 
backscattered  return  from  the  target.  Figure  1  is  an  example  of  a  single  data  record  that  contains 
the  backscatter  signatures  at  each  of  the  four  wavelengths.  Table  1  presents  a  complete  listing  of 
the  times  of  lidar  data  acquisition  and  the  times  of  target  detection  for  16  July  and  19  through 
27  July.  Figure  2  displays  graphically  the  information  listed  in  Table  1.  (All  of  these  exhibits 
are  reproduced  from  the  first  Status  Report.) 

C.  CHARACTERIZATION  OF  RECEIVER 

LOGARITHMIC  AMPLIFIER  RESPONSE 

A  second  task  was  to  characterize  the  response  of  the  logarithmic  amplifiers  that  were  used 
in  the  lidar  system.  This  information  is  necessaiy  for  analyses  of  relative  backscatter  signatures, 
including  proper  removal  of  the  geometric  range  dependence.  During  the  field  tests,  neutral- 
density  filters  of  known  attenuation  (nominally  3,  10,  20,  and  30  dB)  were  inserted  periodically 
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FIGURE  1  EXAMPLE  OF  AN  A-SCOPE  DISPLAY 


in  front  of  the  lidar  receiver(s)  before  recording  data  for  a  number  of  lidar  shots.  For  each  wave¬ 
length,  the  mean  maximum  backscattered  returns  (in  counts)  from  the  target  were  then  plotted  as 
function  of  filter  attenuation  (in  dB).  Linear  least-squares  methodology  was  used  to  derive  the 
optimum  slope  (counts/dB)  of  each  curve.  Figure  3  shows  the  data  and  the  results  of  the  least- 
squares  fit  for  a  calibration  performed  on  21  July.  The  points  represent  the  mean  measured  val¬ 
ues,  and  the  lines  represent  the  regression  fits,  where  saturated  (255  counts)  or  nearly  saturated 
target  return  data  have  been  excluded  from  the  regression  analysis.  Table  2  lists  the  results  of  the 
analyses  of  the  various  calibration  runs.  In  particular,  Table  2  lists  the  mean  and  one  standard 
deviation  counts-to-decibel  conversion  factors  from  series  of  calibration  runs  that  apply  to  par¬ 
ticular  data  files.  The  identification  data  file  numbers  to  which  the  table  values  apply  and  the 
files  containing  the  calibration  runs  used  in  calculating  the  values  are  also  listed.  TTiese  factors 
are  to  be  used  in  converting  the  digitized  signals  in  counts  to  relative  decibel  values. 

D.  TARGET-DERIVED  PATH-LENGTH  TRANSMISSION  VALUES 

Lidar  signals  backscattered  from  the  reflective  target  positioned  928  m  from  the  SRI 
transmitter/receiver  were  analyzed  to  calculate  path-length  transmission  values  at  0.53,  3.8,  and 
10.6  |im  for  selected  time  periods  during  the  FLAPIR  program.  Hard-copy  and  digital  output  of 
these  transmission  histories  were  forwarded  to  Col.  Koenig  for  subsequent  comparison  with 
measurements  from  other  instruments.  Data  were  transmitied  via  electronic  mail.  As  an  exam¬ 
ple,  Figure  4  shows  the  transmission  history  for  the  p>criod  2244  -2322  EDT  on  27  July. 
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Table  1 

TIMES  OF  UDAR  DATA  ACQUISITION  AND  TARGET  DETECTION  FOR  16  JULY  AND  19-27  JULY 
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Table  1  (Continued) 

TIMES  OF  LIDAR  DATA  ACQUISITION  AND  TARGET  DETECTION  FOR  16  JULY  AND  19-27  JULY 
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Table  1  (Continued) 

TIMES  OF  LIDAR  DATA  ACQUISITION  AND  TARGET  DETECTION  FOR  16  JULY  AND  19-27  JULY 
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Table  1  (Concluded) 

TIMES  OF  LIDAR  DATA  ACQUISITION  AND  TARGET  DETECTION  FOR  16  JULY  AND  19-27  JULY 
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Note:  Frequency  of  target  detection  exceeded  95%  for  the  time  periods  listed,  except  for  those  times  mari<ed  by  an  asterisk  (*),  where  the  frequency 

was  between  75%  and  95%,  of  by  double  asterisks  (**),  where  the  frequency  was  less  than  20%. 


DATA  FILE:  CAI>28,OUT 


FIGURE  3  MEAN  LIDAR  TARGET  RETURN  MEASUREMENTS  AND  LEAST-SQUARES  FITS  FOR 
21  JULY  1990 


Table  2 

SRI  FLAPIR  FOUR-WAVELENGTH  LIDAR  CALIBRATION  RESULTS 

(counts/dB) 


Applicable  _ Wavelength  (um) _  Calibration 


Data  Files 

0.53 

1.06 

3.8 

10.6 

Runs  Used 

06-31 

32-46 

47-58 

5.20  ±0.09 
5.09  +  0.02 
4.50  ±0.01 

9.15  ±0.05 
8.51  ±0.04 
6.07  ±0.09 

6.04  ±0.15 
628  ±.009 
426  ±  0.12 

5.56  ±0.09 
5.50  ±0.38 
3.99  ±0.01 

15.22.28,29 

32.41 

54  (first),  58 

V 
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FIGURE  4  TARGET-DERIVED  TRANSMISSION  PERCENTAGES  FOR  2244-2322  EDT,  27  JULY 
1990 


E.  RANGE-RESOLVED  BACKSCATTER  RETURNS 

The  utility  of  the  data  set  for  calculating  range-resolved  backscatter  values  was  investi¬ 
gated  briefly.  Several  examples  of  range-resolved  backscatter  signatures  were  forwarded  to 
Col.  Koenig  in  a  letter  on  19  November  1991.  Some  of  these  are  reproduced  in  this  report. 
Figtire  5(a)  is  an  example  of  range -resolved  lidar  backscatter  signatures  (in  dB)  at  the  four 
wavelengths  for  a  single  shot  along  a  relatively  clear  path  on  21  July.  The  data  have  not  been 
corrected  for  the  range-squared  dependence  of  the  return  signals.  Figure  5(b)  displays  the  same 
data  with  the  geometric  range  dependence  removed.  The  increase  in  signal  with  range  beyond 
about  350  m  at  the  three  longest  wavelengths  indicates  that  the  signal  has  fallen  below  the  level 
of  sensitivity  of  the  system  at  those  wavelengths;  these  data  are  not  valid.  (System  sensitivity 
was  substantially  reduced  to  observe  backscatter  during  fog  conditions,  and  this  explains  why  the 
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LIDAR  SIGNAL 


(b)  AFTER  RANGE-SQUARED  CORRECTION 


lidar  is  able  to  observe  backscatter  only  for  short  ranges  during  clear-air  conditions.)  Figure  6 
shows  analogous  plots  for  a  lidar  record  measured  under  fog-free,  but  somewhat  hazy  conditions 
on  19  July.  The  relative  backscatter  signatures  indicate  that  the  atmospheric  extinction  was  too 
small  to  be  evaluated  at  other  than  the  0.53-nm  wavelength.  Figure  7  shows  the  relative 
backscatter  returns  after  the  range-squared  correction  for  a  relatively  foggy  condition  on  16  July. 
The  data  reflect  large  variations  in  fog  density  along  the  path;  therefore,  it  may  be  difficult  to 
correlate  path-averaged  data  with  in-situ  data.  However,  there  does  appear  to  be  a  high  degree  of 
correlation  among  the  range-resolved  data  at  the  various  wavelengths.  Also,  the  averaging 
effects  of  the  longer  pulse  lengths  at  10.6  |im  and  3.8  pm  are  evident.  Finally,  Figure  7 
illustrates  one  of  the  primary  difficulties  in  deriving  range-resolved  optical  parameters  ftom  the 
data,  namely,  difficulty  in  aligning  the  range  scale  of  the  backscatter  returns  at  the  various  wave¬ 
lengths.  This  problem  is  discussed  in  more  detail  in  Section  V. 
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LIDAR  SIGNAL 


RANGE  —  meters 

(b)  AFTER  RANGE-SQUARED  CORRECTION 


FIGURE  6  LIDAR  RANGE-RESOLVED  RELATIVE  BACKSCATTER  SIGNATURES 
FOR  0219:26  EDT,  19  JUL>^^990 


FIGURE? 


80 


60 


-20 


-40 


0  200  400  600  800  1000 


RANGE  —  meters 


LIDAR  RANGE-RESOLVED  RELATIVE  BACKSCATTER  SIGNATURES 
FOR  0259:41  EDT,  16  JULY  1990.  AFTER  RANGE-SQUARED  CORRECTION 
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5.  DISCUSSION  AND  RECOMMENDATIONS 


A  number  of  issues  were  identified  in  analyzing  the  lidar  data  set  collected  during  FLAPIR. 
These  were  detailed  in  the  November  1991  letter  to  Col.  Koenig.  These  considerations  and 
others  are  summarized  below: 

Variability  in  Background  Atmospheric  Conditions  Between  and  Within 
Tests — Calculation  of  reliable  path  transmission  values  using  the  target  method  requires  an 
accurate  knowledge  of  the  signals  backscattered  from  the  target  under  clear  atmospheric  condi¬ 
tions.  Clear-air  backscatter  returns  from  the  target  were  not  always  readily  available  before  or 
after  the  FLAPIR  runs,  owing  to  changing  atmospheric  conditions.  Hence,  the  clear-air  target 
returns  were  generally  assumed  equal  to  the  0-dB  intercept  value  derived  from  the  least-squares 
fits  of  the  target  return  data  acquired  during  the  calibration  time  periods  when  neutral-density  fil¬ 
ters  were  inserted  into  the  receiver  optical  path.  This  assumption  can  lead  to  large  biases  in  the 
target-derived  transmissions  if  the  assumed  clear-air  target  return  is  incorrect  by  only  a  few 
counts,  and  could  explain  any  gross  discrepancies  (biases)  between  SRI  lidar-derived  path  trans¬ 
missions  and  those  reported  by  other  investigators  using  other  instruments. 

Possible  Variability  of  the  Reflectance  of  the  Target  Due  to  Changing  Weather 
Conditions — ^Target-derived  transmission  values  could  vary  significantly  if  the  reflectance  of 
the  target  exhibited  a  strong  dependence  on  its  state  of  diyness/wemess,  which  is  determined  by 
atmospheric  conditions.  No  quantification  of  this  effectivity  was  attempted  during  FLAPIR. 

Variability  of  Transmitted  Laser  Pulse  Length  Among  Wavelengths  and 
Between  Lidar  Shots — ^This  affects  the  peak  power  transmitted  by  the  laser  and  therefore 
introduces  uncertainties  in  the  target  returns.  Backscatter  smoothing  effects  of  the  longer  pulse 
lengths  of  the  3.8-  and  10.6-^im  lidar  signatures  relative  to  the  0.53-  and  1.06-)im  lidar  signatures 
are  obvious  on  the  data  plots  of  Figures  5  through  7. 

Alignment  of  the  Muitiwavelength  Signals  with  Respect  to  Range — ^This  can  be 
difficult  under  some  atmospheric  conditions  when  calculating  and  comparing  range -resolved 
backscatter  signatures.  The  backscatter  signatures  were  aligned  by  first  identifying  the  location 
of  the  target  returns  within  the  backscatter  signal  array.  This  was  done  by  searching  (separately 
for  the  range-resolved  backscattered  returns  at  each  wavelength)  for  the  maximum  backscattered 
signal  within  a  predefined  range  (i.e.,  in  the  vicinity  of  the  target)  for  each  shot.  The  multiwave- 
length  backscatter  signatures  were  then  aligned  either  by  aligning  the  maximum  target  returns  for 
each  wavelength  or  by  searching  backward  (for  each  wavelength)  toward  the  laser  for  the  point 
where  the  slope  of  the  backscatter  trace  begins  to  increase  significantly  toward  the  maximum  tar¬ 
get  return  and  aligning  the  traces  on  this  point.  However,  neither  of  these  methods  successfully 
aligned  the  backscatter  signatiucs  all  the  time  because  of  varying  laser  pulse  widths.  Under  rela¬ 
tively  clear  air  conditions,  the  maximum  target  retums  are  of  large  amplitude  and  easy  to 
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identify.  However,  the  identification  process  becomes  problematic  as  the  atmospheric 
attenuation  increases  to  the  point  that  the  maximum  target  amplitude  becomes  comparable  to  the 
noise  in  the  signal. 

Large  Dynamic  Range  of  Atmospheric  Backscatter  Vaiues,  and  Limited 
Dynamic  Range  of  Logarithmic  Ampiifiers — ^The  1.06-|jm  wavelength  channel  posed  the 
biggest  problem.  Because  of  its  high  sensitivity  Garge  number  of  ADC  counts  for  1  dB  of  atten¬ 
uation),  it  was  necessary  to  attenuate  the  1,06-pm  laser  signal  by  about  30  dB  during  the  data 
acquisition  to  keep  the  1.06-pm  target  signal  from  saturating  while  at  the  same  time  being  able  to 
retain  sensitivity  to  the  atmospheric  backscatter.  Even  so,  subsequent  calibration  showed  that  the 
target  returns  remained  saturated  until  an  additional  20-dB  optical  filter  was  inserted  in  the  opti¬ 
cal  path.  As  a  result,  at  the  1 .06-)im  wavelength  only,  it  was  not  possible  to  derive  path  trans¬ 
mission  values  using  the  target  method  whenever  the  one-way  transmission  exceeded  about  10%. 

Lack  of  Information  on  Transmitted  Laser  Energy— The  laser  energy  monitors  did 
not  function  properly  during  the  trial;  hence,  no  such  data  are  available.  This  negates  the  ability 
to  normalize  out  shot-to-shot  energy  variations  and  complicates  analyses  of  temporal  backscatter 
fluctuations.  However,  the  laser  energy  variations  are  normally  within  10  percent  and  should  not 
greatly  effect  derived  transmission  histories. 

SRI  recommends  that  these  effects  be  taken  into  account  in  performing  additional  analyses 
of  the  SRI  PTAPIR  lidar  data. 
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